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Abstract. Based on the QCD factorization approach we analyse the branching ratios for the channel 
B pn. From the comparisons with experimental data provided by CLEO, BELLE and BABAR 
we constrain the form factor F^^'^{m^ ) and propose boundaries for this form factor depending on 
the CKM matrix element parameters p and rj. 

o 

(N 

O : 1. NAIVE FACTORIZATION 

m 

^ ■ The investigation of B decays requires a knowledge of both the soft and hard interactions 

^ ^! which control the dynamics of quarks and gluons. Because the energy involved in B 

I I decays covers a large range, from down to Aqcd, it is necessary to describe the 

phenomenon with accuracy. Recently, the BELLE, BABAR, and CLEO facilities have 
^ . been providing more and more data which can be compared with theoretical results and 

hence increase our knowledge in this area. 

In any phenomenological treatment of the weak decays of hadrons, the starting point is 
the weak effective Hamiltonian at low energy [1, 2, 3, 4, 5]. It is obtained by integrating 
5^ I out the heavy fields (e.g. the top quark, W and Z bosons) from the Standard Model 

Lagrangian. It can be written as, 

^^// = ^EVc^MQ(Ai)0/(M) , (1) 

where Gf is the Fermi constant, Vckm is the CKM matrix element, C/(/i) are the Wilson 
coefficients, are the operators entering the Operator Product Expansion and ji 

represents the renormalization scale. In the present case, since we analyse direct CP 
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violation in B decays into pTU, we take into account both tree and penguin operators and 
the effective Hamiltonian is. 
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Vubv:^ (Ci of + C20I) - VtbV^q £ dOi 

i=3 



+ h.c. 



(2) 



where q = d. Consequently, the decay amplitude can be expressed as follows. 
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VtbV;gY,Q{PV\Oi\B) 

i=3 



+ h.c. , (3) 



where {PV\Oi\B) are the hadronic matrix elements, and P{V) indicates a pseudoscalar 
(vector) meson. The matrix elements describe the transition between initial and final 
state at scales lower than jU and include, up to now, the main uncertainties in the 
calculation because it involves the non-perturbative physics. 

The computation of the hadronic matrix elements, {PV\Oi\B), is not trivial and re- 
quires some assumptions. The general method which has been used is the so-called 
"factorization" procedure [6, 7, 8], in which one approximates the matrix element as a 
product of a transition matrix element between a B meson and one final state meson and 
a matrix element which describes the creation of the second meson from the vacuum. 
This can be formulated as. 



{PV\Oi\B) 
or {PV\Oi\B) 



--{V\J2,\0) {P\Ju\B) 

--{P\J4i\0) {V\J3i\B) 



(4) 



where the Jji are the transition currents. This approach is known as naive factorization 
since it factorizes {PV\Oi\B) into a simple product of two quark matrix elements, (see 
Fig. 1). Analytically, Fig. 1 can be written down as. 
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Y,VcKMCi{^){Mi\J2i\0){M2\Ju\B) 

i=l 



■ (5) 



A possible justification for this approximation has been given by Bjorken [9]: the heavy 
quark decays are very energetic, so the quark-antiquark pair in a meson in the final 

state moves very fast away from the localised weak interaction. The hadronization of 
the quark-antiquark pair occurs far away from the remaining quarks. Then, the meson 
can be factorized out and the interaction between the quark pair in the meson and the 
remaining quark is tiny. 

The main uncertainty in this approach is that the final state interactions (FSI) are ne- 
glected. Corrections associated with the factorization hypothesis are parameterized [10, 




FIGURE 1. Naive factorization, where Mi and M2 represent the final meson states. 



11, 12] and hence there maybe large uncertainties [13]. In spite of this, there are indica- 
tions that should give a good estimate of the magnitude of the B decay amplitude in many 
cases [14, 15]. In order to improve the estimate of the hadronic matrix element, we will 
briefly present in Section 2 the formalism of QCD factorization, which is an extension of 
naive factorization. We will see how it is possible to incorporate QCD corrections in or- 
der to include the FSI at the first order in Us into the factorization approach. In Section 3, 
we will list our numerical results for the branching ratios related to the channels B ^ pK 
and B (on. In Section 4, we will constrain the form factor F^^^ and propose bound- 
aries for this form factor depending on the CKM matrix element parameters p and r\ . 
Finally, in the last section we will summarize our analysis and draw some conclusions. 



2. QCD FACTORIZATION 

Factorization in charmless B decays involves three fundamental scales: the weak inter- 
action scale Myv, the b quark mass scale m/;,, and the strong interaction scale Aqcd- It is 
well known that the non-leptonic decay amplitude for B PV is proportional to: 

A{B^PV)ocY,Q{^){PV\0,{^)\B) , (6) 

where we have omitted the CKM factor and Fermi constant for simplicity. The matrix el- 
ements {PV\Oi{il) \B) contain non-perturbative effects which cannot be accurately eval- 
uated. The coefficients C/(/i) include strong interaction effects from the scales Mw down 
to nit, and is under control. The aim is therefore to obtain a good estimate of the matrix 
elements without assuming naive factorization. In QCD factorization (QCDF), assum- 
ing a heavy quark expansion when Aqcd and soft coUinear factorization where 
the particle energies are bigger than the scale Aqcd, the matrix elements {PV\Oi{iJ.)\B) 



can be written as [16]: 



{PV\Oii^)\B) = {P\j,\B){V\j2\0) 



l+Y,rnCi'; + ff{AQCD/mb) 



(7) 



where r„ refers to the radiative corrections in tts and j, are the quark currents. It is 
straightforward to see that if we neglect the corrections at the order a^, we recover the 
conventional naive factorization in the heavy quark limit. We can rewrite the matrix 
elements {PV\Oi{iJ.)\B), at the leading order in AQco/ifib, in the QCDF approach by 
using a partonic language and one has [16, 17, 18, 19, 20, 21]: 

{PV\Oii^) \B) = Ff^^'iml) dx^{x)^{x) +Af ^^(m^.) j'^ dyTl,{y)^p{y) 

+ fd^ fdx ['dyT/\^,x,y)M^)M^)My) , (8) 
Jo Jo Jo 

where 0m (with M = V,P,B) are the leading twist light cone distribution amplitudes 
(LCDA) of valence quark Fock states. The light cone momentum fractions of the con- 
stituent quarks of the vector, pseudoscalar and B mesons are given respectively by 
jc,y, and ^. The form factors for 5 — > P and B V semi-leptonic decays evaluated at 
A:^ = are denoted by F?^^{m^) and A^^^(mp). Eq. (8) can be understood via Fig. 2 
where a graphical representation of the factorization formula is given. The hadronic 




FIGURE 2. Graphical representation of the QCD factorization formula. 

decay amplitude involves both soft and hard contributions. At leading order, all the non- 
perturbative effects are contained in the form factors and the light cone distributions 
amplitudes. Then, non-factorizable interactions are dominated by hard gluon exchanges 
(in the case where the 0{h.QCDl^b) terms are neglected) and can be calculated pertur- 
batively, in order to correct the naive factorization approximation. These hard scattering 
kernels [16, 17, 18, 19, 20, 21, 22], 7^.^ and 7j", are calculable order by order in pertur- 
bation theory. The naive factorization terms are recovered by the leading terms of 7]-^ 
coming from the tree level contributions, whereas vertex corrections and penguin cor- 
rections are included at higher orders of Us in 7^^. The hard interactions (at order ^(oc^)) 
between the spectator quark and the emitted meson, at large gluon momentum, are taken 
into account by T^^ . 



2.1. The QCD coefficients a, 



The coefficients a,- [23, 24], have been calculated at next-to-leading order. They 
contain all the non-factorizable effects at order in a^. In order to clearly separate every 
contribution, the coefficients a; are written as the sum of. 



(9) 



where the first term includes the naive factorization, the vertex and penguin corrections, 
while the second term contains the hard spectator interactions. According to the final 

states, the terms a, have to be expressed for two different cases: case A corresponds 
to the situation where the recoiling meson Mi is a vector and the emitted meson M2 
is a pseudoscalar, and vice-versa for case B. For case A, the coefficients take the 
form [23, 24], 
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where the terms «4,7,fc,<7,fc,fl8,/,fc ^nd afo^^ are. 
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In Eqs. (10) and (11) Vm^V'm represent the vertex corrections, H,H' describe hard 
gluon exchanges between the spectator quark in the B meson and the emitted meson 
(pseudoscalar or vector). 2^ 3' ^mT' ^m^i QCD penguin contributions and 

electroweak penguin contributions, respectively. These quantities are a result of the 
convolution of hard scattering kernels G, with meson distribution amplitudes, 4>. We 
refer the reader to Refs. [16, 17, 18, 19, 20, 21] for more details. Other parameters are 
d = Ci{ii) (in NDR), = a,.(/i) (next to leading order), and Cf = (N^ - l)/1Nc with 
Nc = 3. 



3. NUMERICAL RESULTS 



Assuming that all of the parameters involved in QCD factorization are constrained by 
independent studies where the input parameters related to factorization were fitted, 
we concentrate our efforts on the form factor F^^^ depending on the CKM matrix 
parameters p and rj . In order to reach this aim, we have calculated the branching ratios 
for 5 decays such as 5± p°7t^,B^ p^Tt^,B^ p^7t^,B^ p^Tt^ and5± coTt^ 
where the annihilation and p — 0) mixing contributions were taken into account. All the 
results are shown in Figs. 3, 4 and 5, and the branching ratios are plotted as a function 
of the form factor F^^^ and as a function of the values of p and rj as well. 

By taking into account experimental data from CLEO [25, 26, 27, 28, 29, 30], 
BELLE [31, 32, 33, 34, 35, 36, 37, 38, 39] and BABAR [40, 41, 42, 43, 44, 45, 46, 47], 
and comparing theoretical predictions with experimental results, we expect to obtain a 
constraint on the form factor F^^^ depending on the CKM matrix element parameters 
p and T]. Because of the accuracy of the data, we shall mainly use the CLEO and BELLE 
data for our analysis rather than those from BABAR. We expect that our results should 
depend more on uncertainties coming from the experimental data than those from the 
factorization approach (as opposed to naive factorization) applied to calculate hadronic 
matrix element {p7r\J^\B) since in B decays, l/m^ corrections are very small. 



3.5 



3 - 



2 2.5 



T 

+1 

CO 

CO 



1.5 



0.5 









1 _|_ 1 _|_ L 4 - L 4 - 4 -'- 4. -'- 4 -'^ 








1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

F, (m^ 




FIGURE 3. Branching ratio for p^TT*, for limiting values of the CKM matrix elements (Left 

hand-side). Branching ratio for p^Tr", for limiting values of the CKM matrix elements (Right 

hand-side). Solid line (dotted line) for max (min) CKM matrix elements. Notation: horizontal dotted 
lines: CLEO data; horizontal dashed lines: BABAR data; horizontal dot-dashed lines: BELLE data. 



For the branching ratio B p^n (Fig. 3), we found total consistency between 
the theoretical results and experimental data from CLEO and BELLE. However, these 
results allow us to determine an upper limit (between 0.40 and 0.65) for the value of 
the form factor Ff^^. The weak dependence of the branching ratio on the form factor, 
Ff^^, is related to the tree and penguin amplitudes which are mainly governed by the 
form factor Aq^'' rather than Ff^^. Therefore, this branching ratio cannot be used as 
an efficient way to constrain the form factor Ff^^. Note also that the comparison with 
BABAR data shows agreement between theory and experiment when Ff^^ is bigger 
than 0.5. 

For the branching ratio — > p^;r^ (Fig. 3), CLEO gives only an upper limit for the 
branching ratio whereas BABAR and BELLE do not. Based on this upper limit, the value 
of the form factor Ff^^ must be lower than 0.62. We emphasize that this branching 
ratio is strongly dependent on the form factor Ff^^ and hence provides an efficient 
constraint for the value of Ff^^. For the branching ratio B*^ p^n^ (shown in Fig. 4), 
BELLE, BABAR and CLEO give consistent experimental data. The decay amplitude 
related to this branching ratio is proportional to the form factor Ff^^ and thus allows 
us to constrain the form factor effectively. Requiring agreement between experimental 
values and theoretical results yields a central value for Ff^^ which is about 0.3. Note 
that for these three branching ratios their dependence on the CKM matrix elements p 
and rj is strong. Hence we expect to be able to determine limits for their values when 
more B decay channels are taken into account. 




FIGURE 4. Branching ratio for Z?" p*7r^, for limiting values of the CKM matrix elements (Left 
hand-side). Branching ratio for Z?" p^Tr", for limiting values of the CKM matrix elements (Right hand- 
side). Solid line (dotted line) for max (min) CKM matrix elements. Notation: horizontal dotted lines: 
CLEO data; horizontal dashed lines: BABAR data; horizontal dot-dashed lines: BELLE data. 

For the branching ratio B'^ p°;r° (Fig. 4), BABAR, BELLE and CLEO only give 
an upper limit for the branching ratio. However, the branching ratio does not appear to 
be very sensitive to the CKM matrix elements p and rj . That could help us to obtain an 
upper limit for Ff^^ which is not sensitive to p and T]. We therefore need new data to 
go further in this case. Finally, we focus on the branching ratio B^ (on^, plotted in 
Fig. 5. There is no agreement with the CLEO data for values of the form factor Ff^^ 
lower than 0.25 whereas there is a good agreement with BABAR and BELLE for any 



value of Ff^^. Note that in this case the sensitivity of the branching ratio to the CKM 
matrix elements is bigger than that to the form factor Ff^^ and does not allow us to 
draw any conclusions regarding the value of Ff^^. 

To remove systematic errors in branching ratio data given by the B factories, we 
can look at the ratio of the two following branching ratios: ^{B^ — > p^n^) and 
^(B^ p^TT^) . In Fig. 5 we show the ratio, as a function of the form factor Ff^^. 
All the B factory data are in good agreement with theoretical predictions. The results 
indicate that the ratio is not sensitive to the CKM matrix elements p and rj whereas it 
is very sensitive to the value of Ff^^. Comparison with the data shows that Ff^^ is 
between 0.13 and 0.30 (BELLE), 0.05 and 0.20 (BABAR), and 0.10 and 0.35 (CLEO), 
respectively. Assuming that the value of Ff^^ at = is around 0.30, we have 

^(5° ^ p^n^) ^ 14.2 X lO-*^ and =^(5° p^n^) < 1 x 10"^. 




FIGURE 5. Branching ratio for B con , for limiting values of the CKM matrix elements (Left 
hand-side) . The ratio of two pn branching ratios Umiting values of the CKM matrix elements (Right 
hand-side). Solid line (dotted line) for max (min) CKM matrix elements. Notation: horizontal dotted 
lines: CLEO data; horizontal dashed lines: BABAR data; horizontal dot-dashed lines: BELLE data. 

It has to be pointed out that the annihilation contributions in B decays play an 
important role since they contribute significantly to the magnitude of the amplitude. 
The annihilation diagram contribution to the total decay amplitude strongly modifies 
(in a positive or negative way) the branching ratio B — > p^n^ according to the value 
chosen for the phase ^a- This contribution could be bigger than that of p — (O mixing 
but carries more uncertainties because of its endpoint divergence. We emphasise that 
these two contributions (p — ft) mixing effects and annihilation contributions) are not 
just simple corrections to the total amplitude, but are important in obtaining a correct 
description of B decay amplitude. 



4. FORM FACTOR Ff^^ 

Form factors play a major role in the factorization method (naive or QCDF) since they 
represent the transition between two hadronic states. Their computation is non trivial and 
may carry large uncertainties, depending on models being used. These models include. 



say, QCD sum rules, heavy quark effective theory, lattice QCD and light cone QCD. 
With the available experimental data for the branching ratios, it is now possible for us to 
constrain Ff^^ in a model-independent way in QCDF. 

It has to be noticed that the branching ratios depend on both F^^^ and A'^^^. In Fig. 6 

we show the results regarding the form factor Ff^^{m^p) as a function of Nc^^ , where 
we require that all the branching ratios for 5 decaying into pK and (OK be consistent with 
the experimental data provided by CLEO and BELLE. We have excluded the data from 
BABAR since they are currently not numerous and accurate enough. We have included 
uncertainties from the CKM matrix element parameters p (0.190 < p < 0.268) and r\ 
(0.284 < 7] < 0.366) and we have applied the QCD factorization method where all of 
the final state interaction corrections arising at order Us are incorporated. We emphasize 
that the results are model independent. 

0.7 r 




- 0.4 



FIGURE 6. F^^^ as a function of Nt^^ ■ Plot obtained by comparing theoretical results from QCFD 
with experimental data from BELLE and CLEO for the branching ratios B ^ pn and B — > (on. The plot 
includes the uncertainties from the CKM matrix element parameters p and r\ . 



We found a large common region between BELLE and CLEO for the B decay into 
^^^{in^p) varies between 0.3 and 0.57 and N^J^ can take values 

2- 



pjl. From our analysis, F^ 



from 1.25 to 2.25. Their central values are Ff^^[mp 



0.43 and Nc^^ 



1.75. There- 
reduces one of the main uncertainties in the 



suit obtained for the form factor F^^^(m^ 

factorization process. That obtained for the effective number of colours, Nc^-^ , confirms 
previous analysis where naive factorization was applied for the same decays [11]. 

It is well known that the CKM matrix element parameters p and rj are the main 
"key" to CP violation within the Standard Model. Recall that the weak phase is mainly 
governed by the parameter rj that provides the imaginary part which is absolutely 
necessary to obtain an asymmetry between matter and antimatter. Based on our analysis, 
we are not able to efficiently constrain the CKM matrix parameters p and rj from the 
branching ratios for 5 ^ p;r. In fact, the common region allowed by CLEO and BELLE 
data for branching ratios for 5 ^ p;r does not constrain the parameters p and rj. In 
the analysis we used the values 0.190 < p < 0.268 and 0.284 < 7] < 0.366 [48, 49], 



to which the common region corresponds. However, we can try (as an example) to get 
some constraints on p and rj by only taking into account the central values for the form 
factor Ff^^{mp) and for the effective number of colours Nc'^-^ . According to our work, 
we find the following limits: 0.205 < p < 0.251 and 0.300 < 7] < 0.351. 

5. CONCLUSION 

The calculation of the hadronic matrix elements that appear in the B decay amplitude 
is non trivial. The main difficulty is to express the hadronic matrix elements which 
represent the transition between the meson B and the final state. 

We have investigated the branching ratios for 5 — > p;r,5 — > coTt within the QCDF 
approach. Comparisons were made with experimental results from BABAR, BELLE 
and CLEO. Based on our analysis of branching ratios in B decays, we have constrained 

the form factor, F^^^{m^), and the effective number of colours, Nc^'^ . More accurate 
experimental data regarding branching ratios in B decays will provide more accurate 
results, which will be helpful in gaining further knowledge of direct CP violation in B 
decays. 

This work could be extended to more B decays. It would be very interesting to 
constrain our parameters by investigating channels other than pK for branching ratios 
and asymmetries. By including more channels, we will use more experimental data 
and hence be able to obtain better results for our parameters. In the QCD factorization 
framework, annihilation contributions could be subject to discussions. Clarifying this 
point would be very helpful in obtaining more accurate theoretical predictions. For 
example, it is important to solve the problem related to the end point integral diver- 
gence [16] which is parameterized without any strong physical motivation. Moreover, 
the annihilation contributions have not been included within the QCDF method. To 
obtain a consistent framework, it would be better to find a way to include them within 
QCDF 
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